Confidence intervals for nonhomogeneous branching processes and polymerase chain reactions
نویسنده
چکیده
We extend in two directions our previous results about the sampling and the empirical measures of immortal branching Markov processes. Direct applications to molecular biology are rigorous estimates of the mutation rates of polymerase chain reactions from uniform samples of the population after the reaction. First, we consider nonhomogeneous processes, which are more adapted to real reactions. Second, recalling that the first moment estimator is analytically known only in the infinite population limit, we provide rigorous confidence intervals for this estimator that are valid for any finite population. Our bounds are explicit, nonasymptotic and valid for a wide class of nonhomogeneous branching Markov processes that we describe in detail. In the setting of polymerase chain reactions, our results imply that enlarging the size of the sample becomes useless for surprisingly small sizes. Establishing confidence intervals requires precise estimates of the second moment of random samples. The proof of these estimates is more involved than the proofs that allowed us, in a previous paper, to deal with the first moment. On the other hand, our method uses various, seemingly new, monotonicity properties of the harmonic moments of sums of exchangeable random variables.
منابع مشابه
Exact Statistical Inference for Some Parametric Nonhomogeneous Poisson Processes
Nonhomogeneous Poisson processes (NHPPs) are often used to model recurrent events, and there is thus a need to check model fit for such models. We study the problem of obtaining exact goodness-of-fit tests for certain parametric NHPPs, using a method based on Monte Carlo simulation conditional on sufficient statistics. A closely related way of obtaining exact confidence intervals in parametri...
متن کاملInference for Quantitation Parameters in Polymerase Chain Reactions via Branching Processes with Random Effects
The quantitative polymerase chain reaction (qPCR) is a widely used tool for gene quantitation and has been applied extensively in several scientific areas. The current methods used for analyzing qPCR data fail to account for multiple sources of variability present in the PCR dynamics, leading to biased estimates and incorrect inference. In this paper, we introduce a branching process model with...
متن کاملOn $L_1$-weak ergodicity of nonhomogeneous continuous-time Markov processes
In the present paper we investigate the $L_1$-weak ergodicity of nonhomogeneous continuous-time Markov processes with general state spaces. We provide a necessary and sufficient condition for such processes to satisfy the $L_1$-weak ergodicity. Moreover, we apply the obtained results to establish $L_1$-weak ergodicity of quadratic stochastic processes.
متن کاملHarmonic moments of non homogeneous branching processes
We study the harmonic moments of Galton-Watson processes, possibly non homogeneous, with positive values. Good estimates of these are needed to compute unbiased estimators for non canonical branching Markov processes, which occur, for instance, in the modeling of the polymerase chain reaction. By convexity, the ratio of the harmonic mean to the mean is at most 1. We prove that, for every square...
متن کاملThe Polymerase Chain Reaction and Branching Processes
We construct a mathematical model for the polymerase chain reaction and its mutations using the theory of branching processes. Under this model we study the number of mutations in a randomly chosen sequence after n PCR cycles. A method for estimating the mutation is proposed and the variance of this estimator is studied. We also study the distribution of the Hamming distance between two randoml...
متن کامل